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A Bose-Einstein condensate in a tilted double-well potential under the influence of time-periodic 
potential differences is investigated in the regime where the mean-field (Gross-Pitaevskii) dynamics 
become chaotic. For some parameters near stable regions, even averaging over several condensate 
oscillations does not remove the differences between mean-field and iV-particle results. While intro- 
ducing decoherence via piecewise deterministic processes reduces those differences, they are due to 
the emergence of mesoscopic entangled states in the chaotic regime. 
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Experimentally it is possible to generate precisely con- 
trollable double- well potentials for Bose-Einstein conden- 
sates (BECs) (Ref. [1( and references therein). A future 
goal for this system is the realization of mesoscopic en- 
tanglement When combined with a time-periodic 
potential difference between the two wells, a BEC in a 
double well could also be used to investigate quantum 
chaos @, 0, 0, HI • Another system which is widely used 
to investigate quantum chaos is the quantum delta-kicked 
rotor @, @, [|| . Research on quantum chaos includes top- 
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Often, a mean-field approach within the Gross- 
Pitaevskii equation is applied to describe BECs. Still, 
there are noticeable differences between mean-field dy- 
namics and quantum dynamics: only the latter dis- 
plays the well known collapse and revival phenomenon 
(cf. By time-averaging over several of those oscilla- 

tions, these differences usually disappear. However, pre- 
liminary results [1 51 ] for the periodically driven double- 
well potential indicate that even under time-average, 
mean-field dynamics and quantum dynamics can display 
qualitatively different results in the regime for which the 
mean-field dynamics become chaotic. 

In this Letter, these differences are investigated sys- 
tematically. First, the TV-particle Hamiltonian is in- 
troduced for which the Gross-Pitaevskii equation corre- 
sponds to a driven nonrigid pendulum. If decoherence is 
implemented on the A^-particle level via piecewise deter- 
ministic processes, the quantum dynamics can become 
qualitatively similar to the mean-field dynamics. The 
reason for the remaining differences between both ap- 
proaches is the emergence of mesoscopic entangled states. 

To describe a BEC in a double well with single-particle 
tunneling frequency £1 and pair interaction energy 2hn, 
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where creates (annihilates) a boson in well j; /io mod- 
els the tilt and \i\ is the driving amplitude. Such Hamil- 
tonians have been used for schemes of entanglement gen- 
eration 17, lc| : without the periodic driving, entangle- 
ment has been investigated in BECs [19|,l20|. Other appli- 
cations include high precision measurements, many-body 
quantum coherence [21 , 22 1 and spin systems [23} . 



On the level of the Gross-Pitaevskii equation for the 
above model, a wave function is characterized by the vari- 
ables 6 and 4>, where cos [6»/2] (sin 2 [6>/2]) is the probabil- 
ity of finding the condensate in well 1 (well 2) and exp(z</>) 
is the phase between the two wells. The corresponding 
iV-particle wave- function ("atomic coherent states" [24j ]) 
with all particles in this state reads (in an expansion in 
the Fock-basis In, N — n) with n atoms in well 1): 
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The mean-field dynamics can be mapped to that of a 
nonrigid pendulum [3, l25j |: including periodic driving 
the Hamilton function reads (z = cos (6>/2)-sin 2 (6>/2)): 
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- 2z(^ + ^-sin(H r )) , T = m. (3) 

The experimentally measurable [l[ population imbalance 
z/2 can be used to characterize the mean- field dynam- 
ics. Fig. [T] shows typical Poincare surfaces of section. 
The initial parameters were chosen such that tunneling 
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FIG. 1: Poincare surface of section for the forced nonrigid 
pendulum (the mean-field dynamics (J3| are plotted for vari- 
ous starting points at integer multiples of the oscillation pe- 
riod 2n/u>). Closed loops are characteristic for stable orbits 
whereas irregular dots represent chaotic regions. For a BEC 
in a double-well, the parameters correspond to: 

(a) a tilt of 2/io/fi = 3.0, a driving frequency of u) — 3tt, 
an interaction of Nk/Q = 0.8 and a driving amplitude of 
2/ii/fl = 0.9 (i.e. a one- "photon" -resonance [261] ^ . 

(b) the 3/2- "photon" -resonance with Nk/Q. = 0.1, 2no/Q = 
3.0, uj/Q = 2.08 and 2/xi/fi = 1.8, 

(c) all parameters as in a) except for Nk/Q = 0.3. 



in the driven, tilted double-w ell p otential is enhanced by 
"photon" -assisted tunneling [26| (cf. Ref. j^). If the 
interaction is not too low (Nk/Q ?t 0.4... 0.6), regular 
and chaotic dynamics coexist (Fig. [TJa cf. [HI), for low 
interaction the dynamics are regular (Figs.[TJb and[TJc). 

For the parameters corresponding to the Poincare sur- 
face of section in Fig. [TJa, Fig. [5]a displays the dif- 
ferences between ./V-particle and mean-field dynamics 
by numerically calculating (using the Shampine-Gordon- 
routine [1^]) the time-average of the (experimentally 
measurable |l|) population imbalance (J z )/N (which cor- 
responds to the mean-field z/2): 
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where for (J z )/N — ±0.5 the entire condensate is in the 
left, respectively, right well. Each point represents an ini- 
tial condition ([2]). The differences are small if the mean- 
field dynamics are regular (cf. Fig. [TJa) while they can 
be rather large in the chaotic regime (up to half the the- 
oretical limit, max{|z/2 - (J z )/N\} = 1). Most of the 
deviations between iV-particle dynamics and mean-field 
dynamics in Fig[2ja lie within twice the root-mean-square 
(r.m.s.)-fluctuations of the iV-particle dynamics. How- 
ever, contrary to the preliminary results of Ref. [loT ]. for 
many initial conditions in the (classically) chaotic regime 
the differences can be very small; they are large near the 
boundaries of stable regions. 




FIG. 2: (color online) Quantum dynamics (N = 100) versus 
mean-field dynamics using the parameters of Fig. [TJa. 

(a) The difference of the time-averaged population imbalances 
{Jz)t/N and (z/2)t as a function of 101 2 initial conditions 
(zo, (f>o) in a two-dimensional projection of the resulting three- 
dimensional plot (T = 100/fi). 

(b) The time-averaged root-mean-square (r.m.s.)-fluctuations 
{AJz)t/N of the population imbalance as a function of the 
initial atomic coherent state ([2}. 



In Fig. [2Jb, the time-averaged r.m.s.-fluctuations 
°f {Jz)/N reproduce many features displayed in the 
Poincare section in Fig. [TJa. Note that the values for 
the r.m.s.-fluctuations are well above those expected 
for N = 100 particles in an atomic coherent state, 
sin(0) / '(2y/N) < 0.05, thus clearly indicating that more 
than one atomic coherent state is involved. Bose-Einstein 
condensates of N ss 100 have been realized experimen- 
tally [3fj| . both the validity of the two-mode approxima- 
tion will be better and life-times of mesoscopic entangled 
states will be longer than in larger condensates. How- 
ever, even when the calculation is repeated for N = 1000 
particles, the differences in the chaotic regime remain. 
As the (non-linear) Gross-Pitaevskii equation does not 
allow any superpositions, decoherence should reduce the 
differences between mean-field and quantum dynamics. 

In this Let ter, we use a piecewise deterministic process 
(PDP) (Ref. cf. [Hi]) to model decoherence. To avoid 
to have to introduce decoherence also on the mean-field 
level (the atomic coherent states ^ become orthogonal 
in the limit N — » oo), we use the projection on the atomic 
coherent states [24j : 
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Now, the PDP simplifies to having jumps on one of the 
atomic coherent states ^ after time t with probability 

fVmp = 1 — exp(— at) , a = const. > , (6) 

and Hamiltonian dynamics (TTJ) between jumps. The state 
on which the wave-function is projected is determined by 
the probability distribution 

Pe^dQ^ ?L±l\(^\6,0}\ 2 sm(d)d8d0. (7) 
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FIG. 3: Time-averaged population imbalance (J z )t/N for 
various driving amplitudes /ii in a tilted driven double well 
(2/io/fi = 3.0, u> = 3fi, T = 100/n). The BEC initially is in 
the lower well (zo = 1)- Solid line: (J Z ) T /N for N = 1000 is a 
smooth curve as opposed to the mean-field results depicted in 
the inset, which display chaotic jumps for small changes of the 
driving amplitude. Dots in the main plot: If decoherence is 
included via the PDP-process described around Eq. (|6} with 
on average ~ 5 jumps (a = 1/20) the behavior is closer 
to the mean-field dynamics. Many dots lie in the area defined 
by the curves ({J z )t ± (AJ z )t)/N (dashed lines). 



Figure [3] shows that the PDP can qualitativel y re produce 
the results of the Gross-Pitaevskii equation [33|. With- 
out introducing the decoherence, the qualitative differ- 
ence between mean-field and quantum dynamics are quite 
large; averaging over several PDPs would again result in 
a smooth curve within the error-bars in Fig. [3] While 
BEC-research in quantum chaos often assumes the valid- 
ity of the Gross-Pitaevskii equation 0, H, 0, [H, |34|, at 
least for the model investigated here, only decoherence 
can lead to the chaotic behavior predicted by mean field. 

Furthermore, differences between quantum dynamics 
and mean-field dynamics can also occur in the regular 
regime: Fig. 0] shows that, at least for N — 100, the 
differences can even lie above the result for many initial 
conditions in the chaotic regime (Fig. [2ja). One way to 
reduce the differences is to average over the Husimi dis- 
tribution ([7]) (see Fig. |H cf. Refs. 0, [35| and references 
therein). This decreases the peaks of the differences be- 
tween mean-field and quantum dynamics by a factor of 
2 (in the chaotic regime, the factor can be of the order 
of 5). A perfect agreement cannot be expected as the 
averaged probability distribution on the mean-field level 
is always added whereas in quantum mechanics also de- 
structive interference can occur. 

On the level of quantum dynamics, the differences 
could be due to either a distribution of many atomic co- 
herent states - or maybe even mesoscopic superpositions. 
For our model all mesoscopic quantum superpositions of 
all N particles being either in one quantum state or in 
another can be expressed as a sum of two atomic coherent 
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FIG. 4: (color online) Time-averaged population imbalances 
of quantum dynamics (N = 100) versus mean-field dynamics 
at the 3/2- "photon" -resonance of Fig. [Tjb. 

(a) The difference is plotted as a function of the initial con- 
dition (zq, 4>o) in a two-dimensional projection (T = 100/J1). 

(b) As in a) but the mean-field dynamics are replaced by an 
average over the distribution of initial conditions (0. 



states (see the explanation before Eq. ©): 

\i Jsv )=r,(\6 l ,ct> l )+^\e 2 ,fc)) , 0< 7 <2^ (8) 

If both parts hardly overlap, \{6\, 4>i\&2, <fa)\ "C 1, the 
normalization ij ~ l/y/2 and IV'sp) is a highly entangled 
mesoscopic state (for finite N, the only two orthogonal 
atomic coherent states ([2]) are |0, 4>\) and \tt, 4>%)). In a 
two-dimensional projection (cf. FigOc) such a state is a 
bimodal distribution (for N — ■> oo: two delta-peaks). 

To numerically identify if a given wave- function \ip) is 
in a mesoscopic superposition, we start by searching the 
atomic coherent state |0i, 4>i) for which \{ip\0, 4>) | 2 reaches 
its maximum, m\. Around (0i,</>i) we define the set Rl 
by |(0,</>|0i,0i)| 2 > 10~ 3 (cf. FigElc). As both parts of 
the mesoscopic superposition (jHJ should hardly overlap, 
the second maximum — |(V'|02, </ ) 2)| 2 is searched out- 
side the set Rl . The set R2 is defined analogously to Rl 
by |(0,^|0 2 ,0 2 )| 2 > 10- 3 . The fidelity |(V# sp >l 2 still is 
a function of 7, taking its maximum and excluding large 
overlaps (Rl (1 R2 ^ 0) yields: 
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Yet, this only indicates entanglement if > 0.5. With 



m 2 

Cent = Pfid , C ont < Pfid 

mi 



(10) 



even values of <7 e nt ~ 0.5 can identify mesoscopic super- 
positions (Fig. Oc). In Fig. 0a, the maximum value of 
entanglement (evaluated at r = 5 and 10) is plotted as 
a function of the initial condition (zq, 4>q): within the 
chaotic regime (left), entanglement generation happens 
on faster time-scales than in the regular regime (right); 
for longer time-scales (Fig. [5]b) the entanglement in the 
entire chaotic regime is more pronounced. It reaches par- 
ticularly high values near initial conditions with large 
differences in the time-averaged population imbalances 
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FIG. 5: (color online) Entanglement (|10p for parameters as 
in Fig. [TJa (left column) and as in Fig. [T]c (right column), 
(a), (b): Mesoscopic quantum superpositions were identified 
at times t — 5, 10, 15, . . .; the maximum value of a ent is dis- 
played for 101 2 initial conditions (zo, 4>o) and for (a) short 
times (r = 10) and (b) longer times (r = 100). 
(c) Projection of two characteristic entangled states (with 
maxima mi, 7TI2) on the atomic coherent states (J2j) . Black 
(grey/red) regions: |(6>,0|V>)| 2 > 0.16 (> 0.05). Left: z = 
-0.6, 4> = -2.764601535, r = 80, a cnt ~ 72.3% . Right: 
z = -0.98, 00 = -2.701769682, r = 75, cr cnt ~ 33.5%. 
In the left plot, the large blue/grey circle is a typical set R 
around \0,$) with \{9, 0, 4>)\ 2 > 0.001 (cf. Eq. ©). 



(Fig. [2ja). We obtained qualitatively similar results also 
for other values of driving amplitude and interaction. 

To conclude, generation of mesoscopic entangled states 
can be a signature of quantum chaos for a BEC in a 
periodically driven double well potential. We investi- 
gated the driving near multi- "photon" tunneling reso- 
nances 26] which were recently observed experimentally 
for a BEC in an optical lattice [36| . While decoherence 
can lead to a "chaotic" behavior similar to the predictions 
of the Gross-Pitaevskii equation, the differences between 
quantum dynamics and mean-field dynamics are due to 
the emergence of mesoscopic superpositions. If the mean- 
field dynamics are chaotic, the entanglement generation 
is accelerated and its values are enhanced. 
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